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Abstract. We study inhomogeneities in the distribution of the 
^ excursion sets in the Cosmic Microwave Background (CMB) tem- 

^ perature maps obtained by the three years survey of the Wilkin- 

^ son Microwave Anisotropy Probe (WMAP) . At temperature thresh- 

olds |T| < 90iiK, the distributions of the excursion sets with 
^ over 200 pixels are concentrated in two regions, nearly at the 

O antipodes, with galactic coordinates / = 94°. 7, b = 34°. 4 and 

/ = 279°. 8, b = —29°. 2. The centers of these two regions drift to- 
^ wards the equator when the temperature threshold is increased. 

^ The centers are located close to one of the vectors of £ = 3 mul- 

tipole. The two patterns of the substructures in the distribution 
of the excursion sets are mirrored, with = 0.7 — 1.5. There is 
no obvious origin of this effect in the noise structure of WMAP, 
and there is no evidence for a dependence on the galactic cut. 
Would this effect be cosmological, it could be an indication of an 
anomalously large component of horizon-size density perturba- 
tions, independent of one of the spatial coordinates, and/or of a 
non-trivial slab-like spatial topology of the Universe. 
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1 Introduction 



The power spectrum of CMB has provided essential information 
on the cosmological parameters [HE]. Another important source 
of information are CMB temperature maps which are particularly 
useful for the study of possible non-Gaussianity signals of various 
nature. Among the latter are the alignments of low multipoles, 
claimed as anomalous with respect to Gaussian distribution (e.g. 
PjUjSUSiZljHiS] ) . In the present paper we analyse axial inhomo- 
geneities in the distribution of the excursion sets in the WMAP3 
maps; these can be related to the alignments mentioned above. 
We have estimated the centers of excursion sets at given temper- 
ature thresholds and for a given pixel count. The distribution of 
the centers of the excursion sets revealed inhomogeneities, with 
mirroring features, which we analyse below. Though the Galactic 
and interplanetary contamination can be the major contributor 
here, and hence the non-cosmological nature of the low multi- 
poles, we nevertheless discuss briefly the principal conditions for 
arising of the mirroring from topological properties of the Uni- 
verse. 

2 Excursion sets 

For this analysis we used the 94 GHz (3.2mm) WMAP 3-year 
maps. These maps feature the highest angular resolution, and 
are less polluted than lower frequency channels by synchrotron 
radiation of our Galaxy. The |6| < 20° region has been excluded 
as usual, to minimize effects due to the Galactic disk. 

The algorithm for defining of the center of the pixelized excur- 
sion sets, and of their other parameters, has been described in 
pHlfTT] ; temperature independent ellipticity in the excursion sets 
was found in Boomerang's and WMAP's maps, and had been 
found in the COBE maps as well [12] . The algorithm is based on 
Cartan's proof that maximally compact subgroups of Lie groups 
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are always conjugate |[I3], i.e. due to the conditions of existence 
of the center of mass of the given set of points. For Riemannian 
manifolds the center of mass of the points {xi},i = l,...,n is 
defined as 

X = -(Xi + X2 + ... + Xn), (1) 

n 

or that the mass center is the point where the following function 
is minimized: 

n 

y^Y.d{x,x,f. (2) 
1 

Though for present purposes the temperature weighting is not 
important, this procedure in principle enables such generalization 
assigning a weight Tj to each pixel of coordinate Xi. 

Without loosing generality consider a compact subset A ^ M 
with a mass distribution da on A with normalized mass 1. Car- 
tan's theorem states the existence of the mass center if the func- 
tion 

1 r 

f :me M ^ - d{m, afda (3) 
Z 

is convex, achieves a unique minimum at the mass center of A for 
the mass distribution da. The mass center, then, is the unique 
point of vanishing of the gradient vector field i.e. of the linear 
connection or of the covariant derivative V/(x). 

We used this algorithm to obtain the centers of excursion sets of 
various pixel count and temperature intervals. 

3 The Inhomogeneities 

Our analysis has shown that the distribution of the centers of 
the excursion sets, for pixel counts larger than 100 and 200, re- 
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veals anomalous properties. Namely, starting from the tempera- 
ture thresholds about |r| = 8O11K the centers of the excursion 
sets are located in the Northern and Southern hemispheres, and 
denoted as "A" and "B", respectively; their coordinates for pixel 
counts larger than 200-pixel at |r| = 90 iiK are 

/ = 94°.7, b = 34°.4 (A); I = 279°.8, b = -29°.2 (B), 

so that they are nearly antipodal. 

One can analyze compare the location of A and B to the mul- 
tipoles of CMB anisotropy using the Maxwellian vectors of the 
multipoles jHISlZj. While the coefficients aim are describing the 
representation of the temperature by spherical harmonics 

— Tj, — = L aimYirn{n) (4) 
^ m=-e 



with < a£m >= and for Gaussian functions are defining the 
angular power spectrum 

1 ^ 2 

= o/^ , 1 Wim\ , (5) 



2i+l 



the Maxwellian multipole unit vectors {-O^^'*) | i — !,...,£} and a 
scalar coefficient A^^^ characterize the series expansion of a func- 
tion on a sphere. The algorithm in PJH] was used to compute the 
multipole vectors reducing the problem to finding of the roots of 
a polynomial. 

Figure 1 shows the position of the A and B with respect the 
vectors of multipoles ^ = 1 — 4. It is seen that A and B define a 
direction close to one of the vectors of ^ = 3 multipole, drifting 
towards the equator with the increase of the temperature. 

We also calculated the sum vector of multipoles ^ = 2 — 8 as 
shown in Figure 2 (here the module of each vector was weighted 
by 1/^(^+1)): the open star marks the sum oi ^ = 2 — 8 multipole 
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vectors, the black star denotes those of £ = 1 — 8 i.e. together 
with the contribution of the dipole. We see that centers A,B are 
not near either of them. 



To probe the dependence of the inhomogeneities on low mul- 
tipoles, we removed such components from the WMAP 3 year 
maps using the "anafast" and "synfast" functions from Healpix 
package jH]. 



The increase of the number of pixels in the excursion sets with 
the radius from A and B is shown in Figure 3, upper plot. The 
middle and lower plots in Figure 3 exhibit the same but with 
extracted multipoles £ = 2 and £ = 2,3, respectively. 



4 Mirroring 



Figures 4-6 show the distances of A and B from the dipole, from 
the closest pole of multipole £ = 3 and from the sum of the 
poles of £ = 2 — 8 multipoles, all as a function of the temper- 
ature threshold. The middle and lower panels of Figure 4 show 
the dipole-(A,B) distances from multipoles £ = 2 and £ = 2,3 
respectively. 

The locations of A and B are not close to those of the cold spot 

m- 

The role of scan inhomogeneities and noise is tested by analyz- 
ing in the same way the difference of the maps from independent 
radiometers (A-B). For the studied excursion sets, in the noise 
map (A-B) a signal-to-noise ratio around 4:1 is found; the excur- 
sion sets in this case do not show any of the properties described 
above for the 94 GHz (A+B sum) map. 
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Fig. 1. Drift of the centers A and B with the increase of the temperature threshold 
interval and the multipole vectors i = 1 — 4. 
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Fig. 2. Multipole vectors, ^ = 1 — 8, and the sum vectors, ^ = 1 — 8 (black star) and 
£ = 2 — 8 (open star) . 



6 



Fig. 3. Pixel numbers of the excursion sets vs the radius from A and B (dotted hne) 
at temperature step 1 (upper plot). The middle and lower plots are with £ = 2 
and i = 2,3 multipoles removed, respectively. 

5 Discussion 



The centers of excursion sets have nearly antipodal locations. 
The inhomogeneity centers are close to one of ^ = 3 multipole 
vectors, but not close to the sum vector of multipoles up to f = 8. 
They are close to the ecliptic pole and are nearly orthogonal to 
the CMB dipole apexes; however, they drift upon the increase 
of the temperature threshold interval. The effect weakens, but 
does not disappear completely when the ^ = 2 or ^ = 2, 3 modes 
are extracted. There are visible mirroring properties, even if not 
perfect. 

The association to these inhomogeneities with mirroring features 
either to the ecliptic, or to the dipole or to the 1=3 multipole, i.e. 
to unknown interplay of interplanetary and Galactic foregrounds 
is certainly a central issue. In that case, however, one is led to 
reconsider such non-cosmological contributions in the low multi- 
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Fig. 4. The distance of the dipole apexes from A,B; the rest as is in previous figure. 
The for mirroring is 1.50, 1.48, 4.81, respectively. 



3 Correction WMAP 3yrs W [A+B] Alg:H11101 26.09.07 09:32:40 
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Fig. 5. The distance of ^ = 3 (closest vector) from A,B; mirorring = 0.74. 

poles. For now, this seems one of the strong alternatives. 

What if the described features are nevertheless cosmological? Do 
there exist fundamental mechanisms producing such an approx- 
imate mirror symmetry of CMB temperature fluctuations with 
respect to some fixed plane, apart from a pure chance? The an- 
swer is positive, and the simplest possibility is presented by the 
non-trivial spatial topology of the Universe of the type, i.e. 
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3 Correction 



WMAP 3yrs W [A+B] Alg:H11101 27.09.07 12:30:00 




L(2-8)-CENTER 



Fig. 6. The distance of ^ = 2 — 8 sum vector from A,B; mirroring = 0.78. 

with one of spatial coordinates, say z, being identified: z = z + L. 
This topology may be considered as a limiting case of the Uni- 
verse with compact fiat spatial sections having the topology 
if the identification scales Li,L2 along two other spatial coordi- 
nates are much more than L (the slab topology), see P^BlfTTlfTS] 
for the discussions of quantum creation of the early universe with 
such topology. Another possibility is to have compact positively 
curved spatial sections with the radius of curvature much exceed- 
ing the present cosmological horizon Rhor (so that ^tot — 1 <C 1) 
and with the topology S^/Zn, iV > 1. 

As was shown in jT9], in the case of the spatial topology 
oriented along the z axis, a large-angle pattern of a CMB tem- 
perature anisotropy AT [6, 0) is a sum of two terms. The first of 
them has the exact plane (mirror) symmetry with respect to the 
(x, y) plane, namely 

AT{e,(j)) = AT{7r-e,(j)). (6) 

It originates from the Sachs- Wolfe effect at the large scattering 
surface from density perturbations which does not depend on z. 
The second term represents the remaining part of anisotropy and 
does not have any symmetry at all. However, for uqL of the or- 
der of i?hor or slightly more, where ao = a{to) is the present scale 
factor of a Friedmann-Robertson- Walker cosmological model, the 
latter term should be somehow suppressed since the Sachs- Wolfe 
contribution to it from the last scattering surface comes from 
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perturbations having wave vectors with |k| > 27r/L. That is why 
we can expect a total large-angle pattern of AT to have an ap- 
proximate mirror symmetry in this case which should quickly 



The first direct search for such an effect in [2D] using the COBE 
data with a negative result placed a lower bound of the topolog- 
ical scale gqL > 3h~^ Gpc « 0.3i?hor at the 95% confidence (the 
value of i?hor is given for = 0.3, Qa = 0.7). Using the first- 
year WMAP data and adding a different, "circles on the sky" 
method, this upper bound was raised up to ~ Rhor f3] and then 
even up to « l.SRhor [l21|, see also [221^1121]. 

However, even aoL > 2Rhor (when the circles in the sky method 
does not work at all) does not exclude observability of this effect, 
if the z-independent large-scale part of 3D perturbations has a 
sufficiently large amplitude. It should be emphasized that there is 
no definite prediction of a relative amplitude of this 'zero' mode 
of perturbations compared to a generic mode with kz since 
there is no definite theory of how this topology arises. Moreover, 
the recent re-analysis of this problem for the cubic topology 
[L = Li = L2) (for which there may be no approximate mirror 
symmetry at all) using the three-year WMAP data [25] suggests 
that the lower limit in PT] is too optimistic, so even a^L = 1.2i?hor 
is not excluded in this case (in agreement with the lower limit 
obtained earlier in [26]). Clearly, a lower limit on qqL for the 
topology may not be less than that for the more restrictive cubic 
topology, so gqL > 1.2i?hor seems to be still possible for the 
former topology, too. Note also that a possible mirror symmetry 
of the Universe was proposed in [27j, however, with suppression 
of even multipoles (odd point parity). 

Thus, even though obtaining of secure constraints either on the 
torus topology or the compactification scales needs further data 

* An additional effect worsening tliis symmetry at large angles is a contribution 
from the integrated Sachs- Wolfe effect at small redshifts due to a cosmological con- 
stant. 
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and efforts, tfie large scale non-perfect mirroring in the CMB 
maps revealed above, if cosmological, would already indicate that 
the large-scale ^-independent part of density perturbation inside 
our cosmological horizon is anomalously large. 

We thank Paolo de Bernardis for valuable discussions and help. 
AAS was partially supported by the Research Program "Astron- 
omy" of the Russian Academy of Sciences. He also thanks Prof. 
Remo Rufhni and ICRANet, Pescara for hospitality during the 
final stage of this project. 
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